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Abstract. In this paper we obtain the bispectrum of dark matter density perturbations in 
the frame of covariant cubic Galileon theories. This result is obtained by means of a semi- 
analytic approach to second-order perturbations in Galileon cosmologies, assuming Gaussian 
initial conditions. In particular, we show that, even in the presence of large deviations of the 
linear growth-rate w.r.t. the ACDM one, at the bispectrum level such deviations are reduced 
to a few percent. 

Keywords: Modified gravity, Galileon, Matter Bispectrum, Cosmological Perturbations 



Contents 

1 Introduction 1 

2 Action and field equations 3 

3 Background evolution 3 

4 Cosmological perturbations 5 

5 Linear perturbation theory 5 

5.1 Integral solutions for the growing and the decaying modes of DM perturbations 8 

6 Second-order perturbations 10 

6.1 Solution of the evolution equation for the second-order DM density contrast 12 

7 Dark Matter power spectrum and bispectrum 13 

8 Conclusions 18 
A Gauges in linear approximation 19 
B Source terms for the second-order equations of motion 22 
C Background quantities for the second-order DM kernel 29 



1 Introduction 

Several observations, such as those referring to the magnitude- redshift relation for type-la 
Supernovae (SNIa) [1, 2], Cosmic Microwave Background (CMB) temperature anisotropies 
[3-5] and Baryonic Acoustic Oscillations (BAO) features in galaxy clustering [6, 7], suggest 
that the universe is currently undergoing an accelerated expansion phase, caused by the 
presence of a positive cosmological constant or a more general Dark Energy (DE) component 
or a suitable modified gravity model. Assuming that the matter distribution is dominated by 
Cold Dark Matter (CDM), the simplest model that reproduces this effect and fits present data 
is the so-called ACDM one, based on the existence of a cosmological constant term that fills 
the gap between the matter energy density and the critical one. Even though the presence of 
a cosmological constant term A is fully consistent with General Relativity, its value appears 
too small to be explained by fundamental physics [8] . Consequently, alternative models have 
been explored such as, for instance, quintessence [9, 10], f{R) [11], massive gravity [12] (see 
[13] and references therein). 

In this paper we focus on a modified gravity model obtained in the context of the 
Galileon scalar-tensor theory [14]. This theory is obtained by taking the decoupling limit of 
the Dvali-Gabadadze-Porratti model (DGP) [15]. The Galileon is the most general theory 
containing second-order derivatives in the scalar field with some properties: specifically, in a 
flat space-time it preserves the Galilean shift symmetry {d^ir — >• 9^7r -|- 6^), while in curved 
space-times, adding suitable coupling terms between gravity and the Galileon field [16, 17], 
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the model avoids the Ostrogradski instabihty [18]. In addition, on non-Unear scales the self- 
interactions of the Gahleon field screen the fifth force through the Vainshtein mechanism [19] . 
The essence of this mechanism lies in the non-standard kinetic terms (i.e. Oirdf^Trd^ir), which 
decouple the scalar field from gravity at small scales (r <C ry, where ry is a characteristic 
scale around a matter source, named "Vainshtein radius"). Many literature has recently 
appeared on these models and their generalizations [20-30]. Galileon models have been 
extensively studied at late-times [31-46], during infiation [47-54], and a subclass of these 
models has been already compared by observations [55-57]. 

The main goal of this paper is to study the late-time non-Gaussianities (NG) of the 
matter distribution arising from gravitational instability in the cubic covariant Galileon the- 
ory. It is well known that NG can be classified in primordial and late-time. The primordial 
ones come from non-linearities encoded in the infiationary perturbations [58]; these are im- 
printed in the CMB and in the Large-Scale Structure (LSS) of the universe [59-65], and 
should be constrained by present and future surveys [66, 67]. The late-time non-Gaussianity 
in the LSS is generated classically by gravitational instability, when cosmological perturba- 
tions enter non-linear scales. While a Gaussian universe can be completely described by the 
power-spectrum, the deviations from Gaussianity are encoded in higher-order statistics, such 
as the bispectrum and the trispectrum [68, 69]. 

The interest in studying the dark matter bispectrum in the Galileon model comes from 
the possibility to measure the signature of modifications from standard gravity.^ If this is the 
case, the bispectrum can be used to lift degeneracies among different models giving rise to the 
same observed power spectrum and the same background cosmology. We choose Gaussian 
initial conditions, in order to extract only the late-time non-Gaussianity. In particular, we will 
focus on the dark matter bispectrum calculated at tree- level (second-order perturbations), 
since it gives the leading contribution in the weakly non-linear regime. Even though we 
consider models with important modifications in the background and in the growth rate 
w.r.t. ACDM, we will show that the matter bispectrum deviations that we obtain are less 
than 5%. We think that this suppression is connected with a compensation effect when the 
equation of state is w < —0.8. Our results are obtained by using a semi-analytic technique 
both at first and second-order in perturbations. 

The paper is organized as follows. In Sec. 2 we introduce the action for the Galileon 
model we focus on and the resulting equations of motion. In Sec. 3 we solve the equations of 
motion in a Friedmann-Lemaitre- Robertson- Walker (FLRW) universe. In Sec. 4 we introduce 
the setup used to calculate the perturbed equations. In Sec. 5 we analyze the first-order field 
equations. We also find two integral solutions for the growing and the decaying modes that 
simplify many analytic calculations. In Sec. 6 second-order perturbations are considered. In 
Sec. 7 we compute the power-spectrum and the bispectrum for the DM perturbations. In 
Sec. 8 we draw our conclusions and provide some comments. In Appendix A we derive an 
equation for the first-order DM perturbations in different gauges. In Appendix B we provide 
the source terms of the second-order field equations. In Appendix C we show the coefficients 
of the kernel for the second-order DM fluctuations. 

Throughout the paper we adopt units c = h = G = 1, except where explicitly indi- 
cated; our signature is (— , +, +, +). Greek indices run over {0, 1, 2, 3}, denoting space-time 
coordinates, whereas Latin indices run over {1,2,3}, labelling spatial coordinates. 

^For other works on the dark matter bispectrum within other modified gravity models see [70-73]. 
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2 Action and field equations 

In this paper we study the simplest form of the covariant Gahleon model with a potential, 
which reads [14, 16] 

S=IJ d^^^/=ff [M^iR + 2£. + 2Cm] , (2.1) 



where 



(2.2) 



Here, Cj are dimensionless constants, Mpi is the reduced Planck mass and M is a constant 
with dimensions of mass. Cm is the Lagrangian of a pressure-less perfect fluid (DM) with 
density p and four- velocity u^. 

Varying Eq. (2.1) w.r.t. the metric we obtain the Einstein equations 



G 



M. 



pi 



(2.3) 



where 



rp{m) 
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TT;^ TT-i, Dtt - vr.j^ vt.^Iq, tt'" -|- flf^i, vr''^ vr-afl tt'^ 
The equation of motion for the scalar field vr is 



ci 



The stress-energy tensor continuity equation for the DM component, reads 



{m)^J. 



(2.4) 



(2.5) 
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3 Background evolution 

From Eqs. (2.3) and (2.5) we can study the background evolution in an expanding FLRW 
universe with scale-factor a(r), where r is the conformal time. Let vr = tt{t) be the Galileon 
field at the background level and Pm{T) and Ptt{t) the background matter and the Galileon 
energy density respectively. The (0,0) and the components of the Einstein equations 
read 



a2 



Pm + Pn, (3.1) 



M2 

-^{n^ + 2n') = -p^, (3.2) 

where 7i = a'{T)/a{T) is the Hubble parameter, primes represent derivatives w.r.t. the con- 
formal time T and 

_ CiM^ C2 ,2 , 3C3 ,3 
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P. =— ^ - - IP^" (- - ) . (3-4) 

are the scalar field density and pressure, respectively. The equation of motion for the Galileon, 
Eq. (2.5), becomes 

^ + 3 K + 2^-'] - [2^-" + ^'A = ^ (3-5) 

where, without loss of generality, we have defined = MpxH"^. Here T-Lq is the value of the 
Hubble parameter 'H(t) in a FRLW universe today. 




Figure 1. Background evolution of the Galileon model. In the top panels we plot the deviations of 
the Hubble parameter w.r.t. the Hubble parameter in ACDM and the evolution of the DE density. 
In the bottom panel we plot u'7r(a) = p7r(a)/P7r(a)- The parameter values are: Ci — 1.6, C2 — 0.04, 
C3 = 10~^ (red line); ci — 1.5, C2 = 0.04, C3 = 10"'^ (blue line); ci = 11, C2 — 3.8, C3 = 1 (purple line); 
ci —Q,C2 — 3.6, C3 = 1 (yellow line); ci — 10^^, C2 = 3.3, C3 = 1 (green line); ACDM (black line). 

We have studied the background evolution solving Eqs. (3.2) and (3.5). The initial 
conditions are determined fixing an initial vacuum energy density pniTi) and using the back- 
ground equations in the regime Ptt{t) <^ Pm{T)- The DM and the DE energy densities today 
are 0™(ro) = pm{ro) / {3M^{H'^) = 0.27 and J7^(ro) = p^{to) / {3M^{Hl) = 0.73 respectively 
[4]. We take into account the parameter ci 7^ 0, which is the most general potential term 
preserving the Galilean shift symmetry. It acts as a cosmological constant in the case vr' — >• 0. 

In Fig. 1 we show the evolution of 'H{a), ^-^{0) and the equation of state WT^{a) = 
P-k{o) / P-k{o) for the models we are considering. In the limit ci — )> (green line) we have 
noted that the evolution of the background is C2 and C3 independent. This behavior is 
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expected because if ci (or C3) is absorbed through a redefinition of the Galileon field, C3 (or 
C2) is constrained by the condition r2,r(To) = 0.73. In order to have a free parameter that 
allows to decrease the difference between ACDM and our Galileon models it is crucial to 
impose C3 ~ ci 7^ 0. 

4 Cosmological perturbations 

In this section we give some definitions needed to analyze the evolution of the DM pertur- 
bations on sub- horizon scales [74]. Without choosing any gauge the metric can be written 
as 

ds^ = a(r)2 [_(i + 2^)dT^ + 2ididx'dT + [(1 - 2(l))6ij + Xij] dx'dx^ . (4.1) 

Here the dependence of all the perturbations on both the conformal time r and the spatial co- 
ordinates X is implicit. The symmetric trace- free perturbation Xij and cDj can be decomposed 
as 

tJi =uji + diuj , (4.2) 
Xij =Xij + diXj + djXi + DijX , (4.3) 

where and Xi are transverse vectors (i.e. S^^diUj = 0), Xij is a trace-free transverse sym- 
metric tensor {6^^ Xij = ^^''diXjk = 0) and Dij is a trace- free operator defined by Dij = 
didj — (l/3)(5jjV^. Perturbations of the energy-density and the four-velocity of the DM fluid 
can be written as 

/j(x,T)=p(0)(T)[l + 5(f,r)] , (4.4) 
1 



5^0 + 



V'^[X,T 



a 

We can expand any perturbation up to the desired order in this way 



(4.5) 



^ ^ ^(0) + ^(1) + 1 ^(2) ^ _ i_^(n) (4_g) 

In the following we will drop the suffix "0" . At first order we can safely neglect vector and 
tensor perturbations. In fact the first-order vector perturbations have decreasing amplitudes 
and are not generated by the presence of a scalar field. Moreover, the first-order tensor per- 
turbations give a negligible contribution to second-order perturbations. This result cannot be 
generalized to second-order perturbations, since second-order vector and tensor perturbations 
are generated by products of first-order scalars. 

Perturbing the well-known relation u'^Ufj, = — 1, these useful equations can be obtained 

= -V(^) , (4.7) 

^;(2)0 = _^(2) + 3^(1)2 ^ 2diJ^^d'v^^^ + d^V^^^d'v^^^ . (4.8) 

5 Linear perturbation theory 

At the linear level from Eq. (2.3) we obtain four independent equations, the (0,0), the (0,z), 
the trace and the traceless of parts. These are, respectively, 

2M/ VV« + ^M/V^V^x« - (^6M/^ - ^(i)' - (^2M/7^ - j V^.^^^ 
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= a pmS'^ > + I C27r + 6Afpi n - ^^^^ I ^ Jp^2 2"^' 

M3a2 



, , 9c37r' 7^ \ (1)/ 
+ -^27r + —3-^ vr^ ^ (5.1) 



r — I TT^ ^ -\- 



^(1)" + 2?^x^'^' + - V^x^'^ - - 4^a;(i) + 20(i) - 2V^(i) = . (5.4) 

3 

The equation of motion for the Unear perturbation of the Gahleon field, Eq. (2.5), reads 



C2 



^ M3a2 ^ 1^ ^^'^ + M3a2 + M3a2 J + M3a2 

+ (^-^^" + -JW + ll?3^ j V + ^3^, = . (5.5) 



From the time and the space components of the stress-energy tensor continuity equation, Eq. 
(2.6), we obtain 

5(i)' = 3</<«'-VVi), (5.6) 

^(1)' + ^(1)' + nJ'^ + ^(1) + nv^'^ = . (5.7) 
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There are many ways to decouple these equations . First of all, it is convenient to 
work in Fourier space. From Eqs. (5.4) we can immediately obtain -0*^^^. In the sub-horizon 
(fc^ ^ Ti^) and quasi-static (|0"| < 7i <^ k'^ approximation, the relevant equations 
we need are (5.1), (5.5) and the derivative of (5.6) 



C2 



2c37r" 
APa? 



-/2 



C37r 

M3a2 



M3a2 
1 



7rW + 



C2vr 



2c37r"7r' cavr'^?^ 



M3a2 



,(1) 



6 



(1) 



(5. 



^(1)" + s^'n + e (^0(1) - ^fc'x^'^) - = . (5.10) 

Combining Eqs. (5.8) and (5.9) to eliminate (p^^^ and x^^^ it is straightforward to obtain 




2M3 



Finally, using Eqs. (5.8), (5.10) and (5.11) we are able to single out an equation for the DM 
perturbation 5^^^ 

+ (5.2) 

where G = (SvrMpj) ^ is Newton's constant and we have defined 

2c3 / „ . n, /\ . cs"^^'^ 



a[T} = l Tri-^ + ^^') + . .9 . ■ (5-13) 

C2M3a2 ^ ^ 2c2M6M2ja4 ^ ^ 



The crucial difference between Eq. (5.12) and the one obtained in the ACDM model is that 
the Galileon acts modifying the Newton's constant at late-times. To recover the standard 
Newton's constant it is sufficient to set C3 = 0. On the left hand side of Eq. (5.12) the 
other modification lies inside the friction term {T-L5^^^') due to the evolution of the Hubble 
parameter. As shown in Fig. 1 these differences cannot be neglected and should play an 
important role in the growth of structures. 

Eq. (5.12), which describes the dynamics of DM perturbations on sub-horizon scales, 
together with Eqs. (5.6), (5.7), (5.9) and (5.11) forms our complete set of equations that 
allow to solve the dynamics of the fluctuations at first order. In Appendix A we show how to 
obtain the same result in the Poisson, spatially fiat and synchronous gauges. In particular it 
is important to pay attention doing the sub-horizon approximation in the synchronous gauge, 
due to the residual gauge freedom. 

Eq. (5.12) can be divided in the linear combination of two independent solutions 

j(i)(A?,r) = c+D+(r)5W(fc)+c_D_(r)5(i)(fc), (5.14) 
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where 6^^\k) is the primordial amphtude of the density contrast perturbation. We have also 
added explicitly two integration constants, c+ and c_. D+{t) and D-{t) are the growing 
and the decaying modes and they depend on the coefficients q. In the next subsection we 
will find an integral solution for these modes. 

5.1 Integral solutions for the growing and the decaying modes of DM pertur- 
bations 

To solve Eq. (5.12) it is convenient to redefine 

We can also use the scale factor as the new time variable 

d^6{a) f2 1 dnia)\d6{a) 

+ 7 + ] —T— = Ma)5{a) . (5.16) 



do^ \a Hi^a) da J da 
After that we can perform the change of the variable 



Sia) = uia)J^^ (5.17) 

After a straightforward calculation we shall obtain Eq. (5.12) in its normal form 

d'^u(a) , X , s . X 

- I{a)u{a) = 0, 5.18 

where {—1(a)) is often called the invariant of the equation 

^ , 1 dnia) 1 /dn{a)Y 1 d^n{a) 

J(a) = A(a)+ }' - ——^ + ' . (5.19) 



aTiia) da A'H{a)^ \ da J 2%{a) da^ 
Now, suppose we have to solve 

d^yia) dy{a) dg{a) 

After the substitution 

y(a) = y(a)e+^^<'-'^'^'^("'), (5.21) 
where am is some initial time deep inside the matter dominated era, we obtain 



d'Y{a) 1 
da2 ^2 

We can choose g{a) to be a solution of 



dg{a) 1 2, 



da 2 



9 W 



y(a) = 0. (5.22) 



^-^-\g\a) + 2I{a) = ^, (5.23) 

which is a particular Riccati equation. In this case Eqs. (5.18) and (5.22) become equals. 
Thus, we can relate the solutions of Eq. (5.20) with the ones of Eq. (5.12) through 



y{a) I Uq +1 da'g(a') 

^ ' a \l n{a) 



(5.24) 
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It is straightforward to integrate Eq. (5.20) the first time 



^ + g{a)y{a) = a2, (5.25) 

where 02 is the first integration constant. A second integration is also possible, giving us the 
solutions for y(r) in their integral form 

y{T) = Ki-f\a) + K2-f\a) , (5.26) 

J ami J 

where 

^2(a) = e-/C„'=''''s("'). (5.27) 



From Eq. (5.26) we have two independent solutions of Eq. (5.12) in their integral form 

Di{a) 



7(a) 



a y T-L{a) 
a V n{a) Jam T{a.') 

To determine the growing and the decaying modes it is important to note that there is 
an additional degree of freedom due to the boundary condition in Eq. (5.23). If we want to 
separate them we have to choose carefully the behavior oig{a) at early times. As shown in Fig. 
1, during the matter dominated epoch the contribution of the Galileon field can be neglected, 
leading to an Einstein-de Sitter (EdS) universe. Indeed, during this epoch we expect that 
DEdS^{o) oc a, and DEdS~{o) oc Ti{a)/a oc a~^/^. We can impose Di{a) = D^ds'^iO') = 
obtaining g[a) = — 7/(2a). Taking into account the right coefficients, we can extend this 
result to the general solution, i.e. valid also after the matter-dominated epoch 




a„3/4„y jii^a) \^ 2a„ j,^ ^^{a')) ' ^^'^^^ 

These solutions are important because they are valid in every modified gravity theory 
in which the evolution of first order DM perturbations, Eq. (5.12), is scale-independent. In 
Fig. 2 we show the evolution of our integral solution, Eq. (5.29), vs. the numerical solution 
of Eq. (5.12) for various and arbitrary initial conditions. It is important to note that every 
numerical solution approaches L>+(a), this proves that the first line of Eq. (5.29) is the pure 
growing mode of Eq. (5.12). In Fig. 3 we plot the deviations of the Galileon growth rate, 
/(o) = dlnD/dlna, w.r.t. the growth rate of the ACDM model. For models in which the 
value of C3 is negligible w.r.t. the value ci the deviations are large (up to about 100%), while, 
increasing C3 the deviations decrease reaching ~ 10%. 
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0+ , D 




Figure 2. Evolution of D^{a) (black dashed line), Eq. (5.29), and D{a) (the other lines), solutions of 
Eq. (5.12) with different initial conditions (for a fixed background corresponding to ci — 1.5, C2 — 0.04 
and C3 10~^). 




Figure 3. Growth rate f{a) of the Galileon compared with the growth rate o/ACDM. The values for 
the parameters ci, C2 and C3 are the same as in Fig. 1. 

6 Second-order perturbations 

By perturbing the Einstein and the Gahleon field equations, Eqs. (2.3) and (2.5), at second 
order we can study the dynamics of the DM fluctuations in the weakly non-linear regime. The 
structure of these equations is the same as in the linear case, up to additional source terms 
formed by product of first-order scalar quantities that we will indicate with S'*-"-' (their explicit 
expression in a general gauge can be found in Appendix B). From the Einstein equations we 
obtain, respectively, the (0,0), (0,i) the trace and the traceless part of 
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= a > + I C27r + 6Afpi n - ^^^^ 1 ^ ^p^2 2"^' 

+ aV„V^u,<^> + (c.,' - ^) VV(^) + |lJvV«^>' + S(^) . (6.2) 

^M/V^ _ ^(2) ^ ^(2)' ^ 27^a;(2) _ ^V^x^^^^ + 2Mpi2<^(2)" ^ 4Mpi2^0(2)' 

- (-^.^« - 1^ ) - - - w^) 

+ 5-.mW='-^2^S^+S<^ (6.3) 

V2y2 / (2)" ^ 2Wx<''' + - - 4Woj<'^> + 2,i.(21 - 2*12) j = 251"" . (6.4) 

Eq. (2.5), for the Galileon field fluctuations, becomes 



^ M3a2 ^ 1^ ^ M3a2 + M3a2 ^ ^ M3a2 

/ 2c37r"7r^ 3c37r^\ ^2 (2) ^ C3V'V^a;(2)^ 

+ 1^-^^" + M3^ + j V + ^3,2 = ^ • (6-5) 

The stress-energy tensor continuity equation reads 

5(2)' = 3^(2)' _ v2^(2) + 5(6) ^ (g_g) 
V2 L{2)' + „{2)' + ^^(2) + ^(2) ^ ^^(2)\ ^ ^(7) _ (g^^) 
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In Eqs. (6.2), (6.4) and (6.7) second-order vector and tensor perturbations were present. In 
order to decouple scalar from vector and tensor perturbations we have used the operator di 
in Eqs. (6.2) and (6.7), while we have used didj in Eq. (6.4). Once the equations of motion 
for the scalar perturbations are obtained the steps to obtain the evolution for 5 are the same 
as in the linear case. The result is 

+ M<^V - 4.0 (l - ^-Sgl^) «V™*'=' - S(') . (6.8, 

where 

+ ^(6)' + - S(^) (6.9) 

6.1 Solution of the evolution equation for the second-order DM density contrast 

In this section we study the behavior of Eq. (6.8). It is clear that the homogeneous part of 
this equation is equal to Eq. (5.12). Thus, using Green's method, and Eqs. (5.29), we can 
find an analytical (in its integral form) solution for the evolution of the second-order DM 
density perturbations. Using Eqs. (A. 3), (A. 4), (A. 5), (5.6) and (5.12), in the Poisson gauge 
the Fourier transform of the source term Eq. (6.9) becomes 



2 /4 



5(1) 5(4) 



A;2 



2Mn-2 



'pi 



S'^^\a,k) = J Skid^k25^^\k-ki-k2)lC{aM,k2)5^^\a,ki)5^^\a,k2)- (6.10) 
Here, the symmetrized kernel /C(a, /ci, /C2) reads 

ki-k2){ki^ + k2^) (ki-k2^^ 



/C(a,A;i,A;2) = 71(a) + 72(0) , 2, 2 ^ 73(a) ^ ^ 

ki fc2 ki k2 

a^n^{ki-k2) .2-^2 fa^n^ a^W 

+ 74(a) + 75(a)^ + 76(a) [j^ + ^ 

a^likl + k^ a^n' 

where the background functions 71(a), 72(a) and 73(a) are shown in the next section, while 
the other 7i(a) are listed in Appendix C. Finally, using Green's method with the homogeneous 
solutions, Eq. (5.29), we can find the evolution of the second-order density fluctuations 

5(^)(a, k) = D,ia)6i^Hk) - D,[a) f da'^^Qf^!^ 

Jam a' n^{a')W{a') 



where W is the Wronskian 



5Hr 



+ L>_(a) / ^ / "^v-y- — .g -^2) 



W{a) = D+[a)DJ{a) - D^(a)D+'{a) = "^^2^^ , (6.13) 

am is some initial time deep inside the matter dominated era and 5^'^\k) is the initial second- 
order DM perturbation. It is interesting to see that in this relation there is no explicit 
dependence on the coefficients Cj. 
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7 Dark Matter power spectrum and bispectrum 



To describe the DM distribution of the universe the first statistical interesting quantity is 
the power-spectrum 

(<5(a, h)S{a, fca)) = {27rf d^'^^^h + k2)Pia, h) , (7.1) 

where 5^^\ki + ^2) is the three dimensional Dirac delta function and (. . .) indicates ensemble 
averaging. Note that, under the assumption of spatial isotropy, the power-spectrum depends 
only on the absolute value of ki. By the Wick theorem, for Gaussian distributed fluctuations 
the power-spectrum contains all the information about the DM distribution. The linear 
power-spectrum, calculated using first-order equations, reads 



P{a,k)^\D+{arT'{k)i^—J , (7.2) 

where is the scalar spectral index of primordial fluctuations and T{k) is the transfer func- 
tion (for which we use for simplicity the fit provided in [75]). In the following computations 
we will take = 0.96 [5]. The second statistic of interest is the bispectrum, defined by 

(^5ia,ki)6ia,k2)5{a,k3)^ = {2-nf5^^\ki + k2 + h)B{a,ki,k2,k^), (7.3) 

where the Dirac delta function imposes that only closed triangle configurations are to be 
considered. Since we are interested in studying the contribution generated by gravitational 
instability at late times in the Galileon theory, we impose Gaussian initial conditions. It is 
convenient to use the reduced bispectrum [76], defined by 

n( 1 I. 1 \ - B{a,ki,k2,k3) 

Q{a,ki,k2,k3) = , (7.4 

P{a, ki)P{a, k2) + eye. 

which has the property that it remove most of the scale dependence to lowest-order (tree- 
level) in non-linear perturbation theory. Using the results of the previous sections we can 
write the density contrast perturbation as ^ 

6{a,k) = 6^^\a,k) + ^6^^\a,k) = D+{a)6^^^k) 

+ I d\i j d'q25^^\k-qi - q2)F{a,qi,q2)S^'\a,qi)S^^\a,q2) , (7.5) 



where 



2a''^'H'^{a')W{a')D+'^{a) 



F{a,q,,q2)= I da ^ -,2^2f.,.u,(.,^n ^ ' ^Csuia ,q,,q2) 



7(«) Ja^ \Ja' l\a")) a'^y^W) 



^Notice that in the following we neglect the contribution proportional to the initial second-order DM 
perturbation S^'^''(k) in Eq. (6.12). S^'^\k) contains both a possible primordial NG, and a non-primordial 
contribution, see, e.g. [59, 77]. However the non-primordial term gives a negligible contribution to our final 
results on the scales of the quasi-static regime. 
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The kernel ICsh{cl' ,Qi,Q2) is the leading order of Eq. (6.11) taking into account that we are 
working on scales much smaller than the horizon {ki^ ^ a^H). This kernel can be recast in 
a more convenient form as 



2W{a) (7.7) 



2 



where 



71(a) = / (a) + ^,,0^,0 - 1 + 



72(a) = f\a) + 
73(a) = f{a 



2, . PmO^ cla^7i'^TT''^Pm 



cyn^n'^pi 



Here we introduce h{a) = g{a) + 7/(2a), where g{a) is the solution of Eq. (5.25), to 
parametrize the contribution of the accelerated expansion on the growth rate. Eq. (7.7) is 
one of the main results of our paper. It reduces to the usual form of the Newtonian kernel 
in the limit of an EdS universe [68, 69]. It shows that the different contributions to the 
bispectrum have the same scale dependence as in EdS and ACDM, while they are modulated 
by time dependent coefficients that depend on the particular Galileon model. Looking at 
Eq. (7.7) we can recognize three kind of modifications w.r.t. the ACDM kernel. The first is 
due to the different evolution of the growth rate w.r.t. ACDM and, as stated before, should 
produce deviations in the bispectrum that can reach ~ 100%. The second comes from the 
different evolution of the background, while the third is related to the parameters C2 and C3. 

The reduced bispectrum, Eq. (7.4), assumes the standard form 

. . 2F{a,ki,k2)P{a,ki)P{a,k2) + eye. 

Q{a,ki,k2,k3} = — — — — — — — . (7.9) 

P{a, ki)P{a, k2) + eye. 

The scales at which our approximations can give valid results are IQ-^hMpc"^ < k < 
10~^hMpc~^. The first inequality follows from the sub-horizon approximation, while the 
second excludes the scales at which highly non-linear effects become non- negligible. In Figs. 
4 and 6 we show the angular dependence of the reduced bispectrum for different Galileon 
models, at a = 1 and at a = 0.6 respectively, fixing ki = k2, 6 being the angle between 
ki and k2 {ki ■ ki = kik2Cos6). In Fig. 5 and 7 we show the angular dependence of the 
reduced bispectrum, at a = 1 and at a = 0.6 respectively, fixing ki = eonst. x k2 and 
k2 = 10-3 hMpc-^ 

In Fig. 8 we show the evolution of 

^ _ D+\a') [D.{a)D+{a') - D+{a)D^{a')] JCsnja' ,qi,q2) 
^^"'") = a'^W{a')D^Ha) ^'''"^ 

for an equilateral configuration at a = 1 (left panel) and a = 0.6 (right panel). This 
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ki=k2=0m h Mpc" 



A:,=A:2=0.01 h Mpc" 




0.0 0.2 0.4 0.6 0. 

e/n 

ki=k2 = 0.l h Mpc"' 



e/n 

ki=k2 = 0.l h Mpc"' 




I 




Figure 4. In the left panels we plot the reduced bispectrum for some Galileon models as a ftmction 
of the angle 9, fixing ki — k2 aX a = 1. In the right panels we plot the relative deviations of the 
bispectrum of the Galileon w.r.t. the one of ACDM. The parameter values are: ci = 1.6, C2 = 0.04, 
C3 = 10"'^ (red line); ci — 1.5, C2 — 0.4, C3 = lO^'^ (blue line); ci — 11, C2 = 3.8, C3 — I (purple line); 
ci — 6, C2 — 3.6, C3 — I (yellow line); ci = 10^^, C2 = 3.3, C3 = 1 (green line). 



configuration is useful to understand the behavior of the reduced bispectrum, Eq. 7.9, because 
it is totally independent of the power-spectrum, in fact Q{a,ki,ki,ki) = 2F{a, ki, ki). As 
one can see in the left panel of Fig. 8 the function G{a = l,a') contains a compensation 
effect that reduces the deviations w.r.t. the ACDM model in the bispectrum, as shown in 
Figs. 4, 5, 6 and 7. Let us notice that, for a' < 0.4, the line of every Galileon model we 
consider lies below the ACDM line; viceversa, for a' > 0.4, except for the red and blue lines. 
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Figure 5. The same as Fig. 4 fixing fci = const, x ^2 and k2 — 0.001 hMpc 



for which we find the strongest deviations, the Gahleon hnes he above the ACDM Une (up 
to the present epoch). Consequently, when we integrate G{a = 1, a'), the deviations that we 
have obtained studying the background and the power-spectrum are attenuated considerably. 
Instead, when w > —0.85 - corresponding to the red and blue lines, see Fig. 1 ~ we see a 
minimum below the ACDM around a' ~ 1. This feature decreases the compensation effect 
and produces larger deviations in the dark matter bispectrum. This could be explained by 
the fact that the universe is not accelerating enough today and the evolution of the growth 
rate is strongly modified (see Fig. 3). For these cases the deviations we find in the bispectrum 
are about ~ 5%. Instead, computing B{a, ki, k2, k^) before the acceleration of the universe, 
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0.0 0.2 0.4 0.6 O.S 1.0 0.0 0.2 0.4 0.6 0.8 1.0 




0.0 02 0.4 0.6 O.S 1.0 0.0 0.2 0.4 0.6 0.8 1.0 



e/n e/7T 




Figure 6. In the left panels we plot the reduced bispectrum for some Galileon models as a ftmction 
of the angle 9, fixing ki — k2 aX a — 0.6. In the right panels we plot the relative deviations of the 
bispectrum of the Galileon w.r.t. the one of ACDM. The parameter values are: ci = 1.6, C2 = 0.04, 
C3 = 10"'^ (red line); ci — 1.5, C2 — 0.4, C3 = lO^'^ (blue line); ci — 11, C2 = 3.8, C3 = 1 (purple line); 
ci = 6, C2 = 3.6, C3 = 1 (yellow line); ci = 10^^, C2 = 3.3, C3 = 1 (green line). 

the compensation effect is conserved because the contribution of the Galileon is negligible 
and all models are indistinguishable (see for example the right panel of Fig. 8 and the tiny 
deviations seen in Figs. 6 and 7). 
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0.0 0.2 0.4 0.6 O.S 1.0 0.0 0.2 0.4 0.6 0.8 1.0 




0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 



e/n e/n 




0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 



9/n ein 



Figure 7. The same as Fig. 6 fixing fci = const, x ^2 and k2 — 0.001 hMpc 



8 Conclusions 



In this paper we have presented an analytic expression for the DM modified bispectrum 
in the context of the cubic covariant GaHleon theory. We worked on sub-horizon scales at 
second order in the perturbations, to show the leading contribution in the weakly non-linear 
regime. We have first studied the background with the most general potential that preserves 
the Galilean shift symmetry in a flat space-time. The contribution of C3 is crucial to drive 
the late-time cosmic acceleration, however the deviations w.r.t. the ACDM model are smaller 
if ci ~ C3. At the linear level we have studied the evolution of the DM perturbations finding 
semi-analytical expressions for the growing and the decaying modes. In Fig. 3 we plot the 
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0.2 0.4 0.6 0.8 1.0 0.1 02 0.3 0.4 0.5 0.6 

a' a' 

Figure 8. The integrand of Eq. 7.6 for an equilateral configuration at a = 1 (left panel) and a — 0.6 
(right panel). The models we plot are the same as in Fig. 1, while the black line represents the ACDM 
model. 

deviations of the Galileon growth rate w.r.t. the growth rate of the ACDM model. For models 
in which the value of C3 is negligible w.r.t. the value of ci the deviations are large (until about 
100%), while, increasing C3 the deviations decrease reaching ~ 10%. 

Then, we have extended our analysis to second-order perturbations in order to calculate 
the DM bispectrum. Eq. (7.7), is one of our main results. It shows that the overall k- 
dependence of the bispectrum is the same as in the ACDM model, with time dependent 
coefficients which depend on the particular Galileon model. We noted that, in general, there 
is a compensation effect (see Fig. 8) in the integrand of Eq. (7.7), G{a,a'), that reduces the 
deviations w.r.t. the ACDM model in the bispectrum. This effect is conserved if we compute 
B{a, ki, k2, ks) before the accelerated phase of the universe (see the right panel of Fig. 8), 
because the contribution of the Galileon is negligible and all models are indistinguishable. If 
the bispectrum is evaluated today and the model has w < —0.85, the compensation effect is 
preserved, giving deviations up to ~ 1%. Instead, we noted that this effect is less strong for 
those models which have w > —0.85, allowing for larger deviations in the bispectrum up to 
~ 5%. We argue that the Vainshtein mechanism can be a possible explanation for the overall 
suppression of the deviations w.r.t. the ACDM model in the DM bispectrum and we leave 
for future work further investigation of this aspect. 
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A Gauges in linear approximation 
Poisson Gauge 

This gauge is very useful because in many cases the scalar metric perturbation -0 can be 
interpreted as the Newtonian potential. It can be obtained suppressing the off-diagonal 
terms of the metric 
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= 

a;(^)=0. (A.l) 
From Eq. (5.4), we obtain the standard result 

V'(^) = (t)^^^ , (A.2) 

while Eq. (5.7) reads 

^(1) = 0. (A.3) 
In sub-horizon approximation the field equation for the galileon, Eq. (5.5) reads 

2C37r" ^C3^7r'\ . (l) Cg^'^ 



C2 



A;Vi)-^A:Vi)=0. (A.4) 



M3a2 M3a2 J u'^a? 
Using also the time-time component of the Einstein equations, Eq. (5.1), we obtain 

fi + t". . - ^ - ^) 

Substituting Eqs. (A.3) and then (A.2) into the derivative of Eq. (5.6), in sub- horizon 
approximation we obtain 

5{i)" + ^<5(i)' = _fc2^(i). (A.6) 
Using Eq. (A. 5) to eliminate the metric perturbation (j) in Eq. (A.6), the result is 

(2/4 \ 
1 ^2 .(1) (j^ J^ 

This equations studies the dynamics of the DM perturbation 5^^\ and it is the same 
equation obtained without choosing a gauge, Eq. (5.12). 

Spatially Flat Gauge 

The spatially flat gauge can be obtained by considering the spatial scalar fluctuations equal 
to zero 



<A(i) = 

x(i) = 0. (A.8) 

In this gauge Eq. (5.7) remains the same, while from Eq. (5.4) and its derivative we can 
solve for ip^"^^ 

^(1) = _ 2Uuj^^^ . (A.9) 
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To separate the gahleon perturbation we use Eq. (5.5) in sub-horizon approximation 



Using the last equation in Eq. (5.1), after a sub-horizon approximation, to eliminate 
the galileon field vr^^^ we obtain 



C2 



2c37r" 2c3?^7r'\ 2 



M3a2 M3a2 j " ^™ 



- (^2c2Afpi^ + - ^^^-^ _^^_^JA;^ . (A.ll) 

Finally, to find the dynamics of 6^^^ we have to substitute uj^^'^ from the last equation 
in the derivative of Eq. (5.6). It is straightforward to show that also in this gauge the result 
is identical w.r.t. Eq. (5.12). 

Synchronous Gauge 

The synchronous gauge is a gauge that, at first order, leaves only the spatial scalar pertur- 
bations 

V'(^) = 

J^'>=0. (A.12) 

It is slightly different from the other gauges described, because it has a residual gauge 
freedom. From Eq. (5.7) we find that the velocity v^^^ must satisfy 

+ = . (A.13) 

One can fix the residual gauge freedom imposing the additional condition u^^^ = 0. 
However we do not need to fix it to decouple on sub-horizon scales the DM density fluctuation 
6^^\ Taking the difference between Eq. (5.3) and Eq. (5.1), and performing a sub-horizon 
approximation the result is 

6M,V(^)" + (6M,^,^ + ^) = aW6^'^ + . (A.14) 

In this gauge Eq. (5.5) reads 

3c37r^ (1)// /6c37rV' gc^nn'^ „ A ,(1)/ 

- - f c. - ^ - /tV(i) (A 15) 

Combining Eqs. (A.14) and (A.15) to eliminate the galileon field vr^^^ we obtain 



+ W A« 
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2c3{7r" + 'Htt')\ 2 , ^ 

C2 + ^-j^ j a^p,nS^^> . (A.16) 

It is now straightforward to use this equation, Eq. (A. 13) and the derivative of Eq. (5.6), 

6W" + ^^(1)' = 3 ( 0(1)" + ncP^iY) , (A.17) 



to obtain Eq. (5.12). 

B Source terms for the second-order equations of motion 

In the foUowing we give the exphcit expression in a general gauge of the source terms found 
in Sec. 6. They reads 

12 o y 

+ 4Mpi2</,(i)'v2..(i) - -Mpi^V^x^i^'v^L^^i) - -Mpi^T^V^x^i^V^^^i) +Mpi2v2a.(i)' 



+ ^Mpi2v2;t(^)v20(i) - 24MpiV(')'^0('^ + ^Mpi^V^V^x^^V^'^ " SM^i^nV^J'W^'^ 
3 3 

+ IGMpi^vV^^V^^^ - 4c27r'7r(^)V^^^ + 24MpiV(i)'?^V^^^ + SMpi^^V^w^^V^^^ + 4c27r'V^^^^ 

^94A^ 2..2,(i)2 , 18c3vr-^vr(i)V(i)^ 18C3^V(^)^'?^ C3^-^V^X(^^V^X<^^ , 4c3vrV(i)VV(i) 
+ 24Mpi « ^ + ^^3^2 M3a2 3M3a2 + M3a2 

2c37r^'V^X^'^V'vr(^) 6c37r^'7r(^) V^a;(i) 2c37r^'V^x^^V'c^^'^ 
^ 3M3a2 ^ M3a2 + 3M3a2 

^ 12c3V'(/>(^)V(^) ^ 4c3V'V^7r(i)0(i) ^ 4c37r'2v2a;(i)0(i) 24c37r'=^(/>{i)V{i) 

+ 



M3a2 ^ M3a2 ^ M3a2 M3a2 
72c3A(i)^1^V>(i) icsTT'^V^TrW^pm Scstt'^V^co^^^^^^ 
M^a? M^a? M^a? 



fiA^ 2^2;, (1)^^ , i2c3^^^i^9,^(i)a-^(i) 2cy%^^m^ 

-6M,,nd.ui^dJ^+ M3a2 AW 

+ 4Mpi27^9,t^(i)aV(^) - ^"'''^'^'^^2^^'^^'^ + 6Mpi25,0(i)9V(^) + 4Mpi27^a.a;(i)a>(i) 
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+ -^M^^d^d^d'x^^^d'^djda^^^ (B.i) 
o io o y 

o o 

3 

+ SMpi^V^^d) V^^^ + ^Mpi^V^V^x^^^V^'^ - 8Mpi2v2,/>(i)V(i) - ^Mpi^V^V^x^^^V^^^ 
3 3 

. /w2(i),m 2.,v.2,a) (1) Wtt^W^)' 6c37r'V)Wi) 

12c37rV(i)^^V27r(i) 2c37rW(i)' 203^ V27r(i)v2a;(^) 6c37r^\(^) V2V>(^) 
M3a2 ^ M3a2 + M3a2 M3a2 

_ 8c37r^'V27r(i)V(i) 24c37r^'HV27r(i)V.(i) 12c37r^'vV(^)V>W _ 2 2 (1)^2 (1) 2 

- 2p.V2.(^)^(^)a2 + 4c3ya^(^)-av(^)- ^ _ 12c3.-H^.a),v(^)- 

+ '''''^'m3T'^^'^^ + 16Mp,2a,0(^)aV(^)' - 12Mp,2a,^(^)aV(^)' - ^Mpi^S^V^x^^^^^V^x^^)' 

+ 4Mpi2H5,a;(i)a'V2a;(i) - ^^^^^'^^3^^'"^^'^ - 4M^,' d,J'^ d^V' ej^^'^ 

- 16Mpi2H5,,/;(^)aV(^) + ^^^^^^'^f^^^^'^^'^ + 2pma'd,S^'^d'vW - 4p^a'd,<p(W^ 

- 2p^a^diil}^^'^d^v^^^ - M^i^djd^cp^'^'^d^dix^^y - M^^^djd'i^^^^d^ dix^^'^' + 2M^i^djd\^^'^d^di^'^^y 
+ ^Mpi29,.5'x(')5^9.V2x(^)' - ^Mpi^S.a^x^^^Vs.V^x^') - \M^?d,d^u^'^d^d,V\^^ 

+ 2pma'd,d\^'^d^dM'^ + lM,,'d,d^d\Wd%d,x^'y - ^^"'""^C'^'"^'^' (B.2) 
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(3) _ 12C3^(^)^V(^)^ 6c,n"^'-'y' lSc,7r'\W'^,m' 
IScstt'ttW^H 24c37r(^)'V^^V>(^) 48c3^'VV(i)V(i) 



Af3a4 j^s^i ]vpa 



4 

//„/2„/,(l)2 



M3a4 M3a4 M3a4 

^ M3a4 a2 a2 + ^ 12a2 

2Mpi^V^X^'^'V'x^'^ 4Mpi^HV^x(^)V^x(^) 4Mpi^V^a;(i) V^x^'^ 

Sa^ Sa^ 
Mpi^V^V=^X^^^V'x^^^ 4MpiV(i)V2a;(^) 4Mpi^V>(^) V2a;(i) Mpi^V^x^^) V^a;(i) 

8Mpi^?^V=^X^'^V'^^'^ Mpi=^V2a;(i)' 2Mpi^V2x(^)V^0(^) 4Mpi^V2x^'^ W^^'^ 
_ 24Mp,2</,(i)'V(i) _ 48Mpi2</,(i)'H0(i) _ SMpi^V^a;^^^^') ^ SMpi^V^V^x^^)^^^) 



a 



2 



16Mpi^1^V^a;(i)0(i) 16Mpi^V^(/.(^)(^(i) 8Mpi^V>(^)0(i) 24MpiV(^)>(i) 
12c27rV(i)V(i) 48MpiV(i)^?^V;(^) 48MpiV(i)^7^V;(i) 8Mpi^V^a;(i)V(i) 
16Mpi^HV^a;(^)V>(^) 8Mpi^VV(^)i/>(i) _ 12c27r^V(^)' ASM^^'^H' 4^^^^^ 



24Mpi^HV(^)' 8c3Va,^(i)av(^)' ioc37r^'a^c^(^)av(^)' 

^ a2 ^ M3a4 + M3a4 

6c37r^^Oi7r(^)a^cj(^)' 6c37r^^0ia;(^)a^a;(^)' 12Mpi^-^aia;(^)g^a;(^)' 

^ M3a4 ^ M3a4 

4Mpi2ai0(^)a^w(i)' 12Mpi29iw(i)aV^^^' 2Mpi25iu;(i)a^V2x(^'' 

~^ 2 ' 2 ^ 2 

8Mpi25ia;(i)'5^V2x^^^ _ BMpi^^iV^x^^^^'V^x^^^ 16Mpi2H5ia;(i)a^V2x(^) 
C20i7r(^)a'7r(^) 12c37rVa^7r(^a^a;(^) _ 12c37r^^Hgj7r(^)av(^) 



o2 M3a4 M3a4 

^...„ ....... ...... 



3c27r'^5ia;(^)9^a;(^) 12Mpi2-^'aia;(i)9^a;(i) GMpi'^'H^diUJ^^^d'u^ 



a? c? a? 

SM^.rnO.J^'^O'cM^ 6M^i-(),(:M^0'o^'^'> 2c3V-(),^^i^(>t'(i) 
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_ 4:C37r'^diUj^^^d'ij;^^^ ^M^x^UdiU^^^d'tp^^^ AM^i^di(t)^^'^d'ip^^^ 

Mpi'^djd'u^^^dWiX^^y 4:Mpi'^djd'x^^'>dWiu;(^y _ M^i^djd'x^^^d^ diV'^x^^^ 
a? Sa^ 

~^ 2 2 2 

^ AMpi'^djd'x^'^d^dii^i^^ ^ Mpi-'OkdW^X^^^d^diX^^^ ^. 

H ^ 1 [o-o) 



_ 1^2 (iyy2y2 (1)' _ 2 ,(l)"^2y2 (1) _ uW-Uy^y^ (1) _ 1^2 (1)"^2^2 (1) 

3 6 

3 3 36 

+ Iv^X^'^ V^V^a;^!) + V20(i) Vx^^) + 27^VV^'^ V^x^^^ - iv^V^x^'^'V^x^'^ 
6 6 

- l^V^V^x^i^V^x^i^ + Iv^V^u^^y^V) - J-V^V^V^x^^V^x^i^ + ^7{V2v2a;(i)V 
3 3 18 3 

\2 



3 ■ ■ r ■ ^ ■ 3 ■ ■ " ■ ^ M -2 



(1) 



3 ^ 3 ^ Mpi^ 

3 

- 5V2</.(i)' - v2x(i''VV^'^ - 2?^V2x(i) VV^'^ + 2V2a;(i)'v2V;(i) + - V^V^x^'' V^^'^ 

3 

- 2HV2a;(i)vV^'^ + W^^'^V^^'^ - V^V'^^)' - 2V^V^x^^^" <t>^^^ - mV^V\^^y (t>^^^ 

3 

+ 2V2v2x^^^'V^^^ + 4-HV^V2x^^^V^^^ - ASj'^V'^uj^^y iIj^^^i - SV.V^V^J^^'ip^'^^ - 4V2V^V^^V^^^ 
_ 6c3^Vyi)V7rW 3c3^^V^7r(i)' 6c37r^?^VV(^)' Scg^^'vV^^) V^a;(^) 

9c3V'?^V^7r(i)V^a;(i) 6c37r^'V^7r(i)V^V;(i) 3c37r^V^a;(^)V>(^) 
Mp\ Mpi 6 

- 4dicf>^'^d'v\^'y' + Adii,^'^d'v\^'y' - 2d,ci)^^yd'v\^^y + 2d,'ii^^^y d'v\^^y 
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6 3 3 

M"*Mpi a"' 

3 18 3 



2 



M^MpX^a? Mpi^ M^Mpx^a 

M^Mpi^a^ Mpi^ ' 

6 

M3Mpi2a2 ^ ^ ^ ^ Mpi^ 



- 36>ja^x^^^^'5i'^^^^" - QHdjd'(f)^^^d^dix^^y + endjd'^^^^d^dix^^y 

+ 86>jaV^^^^'5iW^^^' - 6djd'i;^^^dWM^y - ^djd'x'-^y&di^^y 



mdjd\^^^&di(t)^^y - 3djd'J^^d^di<j)^^y + Idjd'x^^yd^dii^ 



2 



- Sdjd'J^^d^di^^^y + ^djd'x^^^&'diV'^x^^^" + djd'x^^ya^diV^^^y 

+ Hdjd'x'^^'>d^diV\^^y - djd'J'^^d^diV\^^y - 2djd\^^'^d^diV'^J^y 
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3 3 3 



+ 



+ 



PmO^djd'-V^^^dWiV^^^ 1q o7oi (l)afcfl (1)" , Ifl «7ai (l)'afca a (1)' 

2 + n^kd^d 'd djdix^'' + -dkd^d'x^ ' d^djdiX^'' 

+ 'Hdkd^d'x^^^d''djdiX^^^' - dkd^d'u^'^^d^djdix^^^' - ^dk& d\^^^d^djdiuj^^^' 
- ^dkd^d'x^^^d'^djdiV^x^^^ - mdkd^d'x^^^d'^djdiJ^^ + dkd^d'J^^d'^djdiJ^^ 

hkd^d'x^'-^d^djdici)^^^ - ldkd^d'x^^^d%di^^^^ - jdid'^d^d'x^^^d^dkdjdix^^^ (b.4) 



+ 2 



^ ^ i»C2« ^3^2 ^3^2 j + TT ^4C2 ^3^3 j 

+ 7r iZC2 ^3^2 ^ ^"3^2 ^^3^2 M3a2 ^ 

+ 7r V a; ^ZC2 ^3^2 M^a^ ' ^ ^ ^ ^ ' M3a2 M3a2 ^ 

+ ^2 (i)^2^(i).2 _ 4c3^^^ _ 4c37r^ ^2^(1) , (i).8c37r^ 8c37r^ 

+ ^(^)^-8c2." - 16C2.'H + + l^^gj^) + ^(^)V^^)(12C2.' 

_ 2403^^- _ 36c37r^ v2a;(^)cA(^)(4c.7r' - - ^^^l 

M3a2 M3a2 > ^ ^^^27r ^3^2 M3a2 ' 

+ V\«VVl)(-C27r' - - 2C3^''^ . 2 (l)^^2 (1), 1 / 
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M3a2 ^ M3a2 ^ M3a2 ^ M3a2 

C37r^V2^(l)^2 4^^y^(l)V2^(l)^ C3VV^X^'^V'X^'^ 

^ 3M3a2 M3a2 ^ 3M3a2 

2c37r^'v"a;(^)V"x(^) 4c37r(^) ^2^(1) 8^^^>^{iy^2^{i) 

3M3a2 M3a2 + M3a2 

4c37rV(i)V27r(i) 4c37r(^^1^V27r(^) 2c37rV2^(i) V2^(i) 2c3VV(^)' 
^ Af3a2 M3a2 3M3a2 ^ M3a2 

4c37r(i)\^V2a;(i) 8c37r^'./>(^) V^a;(^ GcgTr^V^^ V^a;(i) 
M3a2 ^ M3a2 ^ M3a2 

_ 2c37r;VV^)Vyi) 4c37rV2;,(i)v2^(i) 2c3V'V^a;(^)' 
3M3a2 ^ M3a2 ^ M3a2 

4c37r^7r(^)V2V;(i) 2c37r^V^x^^VV('^ 12c30(^)"7r^'(/>(^) 
M3a2 3M3a2 M3a2 

4c3VV^a;(^)V(^) 4c37r^V^V^(^)(/.(^) 24c30(^)\^V(^) 
M3a2 M3a2 ^ M3a2 

8c3V'V^a;(^)V(^) 8c37r^'vV(^)^(^) icadjTrm'd^Trm' 
^ M3a2 ^ M3a2 ^ M3a2 

_ 8c3Ha..(^)av(^)^ 4c.5-^(^)5V(^)' - 24c3V^a.c.(i)gv(^)- 

M3a2 +^C2a,W OTT ^3^2 

8c3Va,V.(^)av(i)^ , ^ . (1)., 12c37r^-^0.7r(^)a-cc;(i)^ 
M3a2 M3a2 

(1)^^ (1)' 18c37r''Haia;(^)5^a;(^)' 2c37r''5i0«a^a;«' 

+ 2c,n dM )a 0.0 ^^3-^ + 

8c3Tr'diTT^^^d'(l)(^y 8c3Tr'^diJ^')d'(j)^^y ic^TT'dinWd'y^x^^^' 
^ M%2 + Jp^ + 

4c3V=^a,a;(^)a-V2x(^)' 4c37r^^a,a;(^)^0-V2x(^) 4 , 2 (1) 

^ 3M3a2 + 3M3a2 3C27ra,a; avx 

^ 3M3^^2 + Jp^ 

3M3o2 M3a2 
12c37r'^'a.7r(i)a^a;(i) 12c37r"^a.^(^)5^a;(i) , ^ n) 
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+ AW 2c,vr )a + 

+ ^^3^^2 2c,5,a )5^9.vr( + 



+ W M3a^ 2c,nd,dx^^d^dM^ 



4:C3iT"iT'djd'x^^^d^diOj^^'> Qc^'K''^ndjd'x^^'^d^diOJ^^'i 2c3^^'^^j^^x'^^'>^^^^iJ'^^'> 
^ ^ M3a2 + M3a2 ^ 

3 

- 2diu^^^d'v^'^y - idM^^d'v^^y - 2di5^^^d'v^^^ - 2ndiJ^^d'v^^^ + 6di(p^^^d'v^^^ 

- 2HdiV^^^d'v^^^ + djd\^^^&dix^^^' - ^diV^^^ d' u^^'^' - Adiip^^^d'v^^^ (B.6) 



3 

+ 2 v2^{i)'v2^(i) + -?^v2x(i)v\(i) - 2V2a.(i)v2r;(i) - 2V\(i)' - 2V2a.(i)'5(i) - 2V\(i)'<^(i) 
3 3 

- 2?^v2^(i)5(i) - 2V2V'(i),5(i) - 2?^v2^;(i)5(i) + W^t^^^^V^'^ + mV^^^^'^^^ + AV^J^^' 4,^^^ 

+ 4a,^(i)a*u;(i)' + 6a,..(i)av(')' + imv^^^d'<i^^^^' + 2d,J'^d'i^^'y - 2d,6^'^d'v^'y + 2v\^'yijW 

3 3 

- ^mM'^^d'V^x^^^ - 2diu;^^^d'V^v^^^ - idM^^d'V'^v^^^ - 2'Hdi5^^^d'u^^^ - 2di5^'^^d'ij^'^^ 
3 

+ mdiJ^^d'ip^^^ + 49^1/^(^)9^^^^ - 2ndiS^^^d'v'-^^ + 4-?^9i(A(^)9'w(i) + 2'Hdii^^^^d'v^^'^ 

- 2djd'v^^^d^diX^^y - 2djd'x^^^d^div^^y - 2'Hdjd'x^^^d^d^v^^^ (B.7) 

C Background quantities for the second-order DM kernel 

In the following we give the explicit expression for the background functions 7i(a) found in 
the kernel (6.11). They reads 



2Mpi2 64c23Mi8MpiiOa3?^2Qio 8c23Mi5Mpi8a4?^2Q8 
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^ TcsVmV''^ Scs^p^'^'V aScsVmV^ 



+ 



l6c2^M^^Mpfa^na^ ^c-i^ M^^ M^f a^Ha^ 2c2^ M^"^ M^x^ a'^ 
C2^M^'^Mp{^a^a^ 2c2^M^^Mpi^a^a'^ 2c2^M^^Mpi^a^n^a^ 

4„ 2^/8 r,f2nj2 



2 



32c22Mi2Mpi«a3H2a6 16c2'^M^^Mpfa'^n^a^ a 



C3^PmW^f Cs^PmTT'^f W PmT^'^ f 

C2^M^^Mpi^a'^Ha<^ C2M^Mpi^aHa<^ C2M^Mpi^aa'^ 



+ 



Wpm^T^'' ^WpJt^''' C3'prna'7r''f 

C22M9Mpi^a3^a4 8c2^ Mpi'^ a^Tia"^ C2M^Mp\'^a'^na* 



+ 77 , , . fi 9o ,9 9 ~ 77 fi 0,9 9 



2c2M^Mpfa^V?o? l%C2M^Mpx^a'^V?a'^ mc2M^Mpi^aV?a'^ 



75(r) = — + 



2Mpi2 Mc2^M^^Mpx^^a^V?a^^ %C2^M^^Mpfa'^V?a^ 

33C35p^2y9 3c3V^^yV 15C3Vm'7r'^ 



16c23Mi5Mpi«a3^Q8 8c23Mi5Mpi«a3Ha8 32c22Mi2Mpi«a3^2(ji 



6 



+ 



4c22Mi2Mpi«a2^2(j6 4c2M6Mpi^aa4 ' 4.C2^M^Mpi^am'^a'^ 

33C3=^P^V'^ 3C3=^p^2^'V 15C3VmV'^ 



8c22M9Mpi6a27^a4 4.02'^ M^i^ ama"^ \Qc2M^Mpi^a'^V?o? 



I5p„j l^Pmf 9C2PmVr'^ , 9C2/9mVr'^ , 9c2apm7r'^ 



^ 2Mpi2 4Mpi2 2Mpi4^2 ^ A^Mpx^cM? 4Mpi^^2 



+ 



2Mpi^V? ' 128c23Mi8Mpii°a3-^2Qio l6c23Mi5MpiSa4^2a8 

ni 4 /10 6 3 /lO 1 c'T 5 2/9 



16c2Mi2Mpi«a^2(j8 l6c24Mi8Mpi^'Ja4^2a8 32c23Mi5Mpi«a3^a8 

.5^ 2^/9-. r-A^ 2^/6 2^/6, 



, 9C3V^^7r^7 , 5C3>^V'° 5C3Vm^7r'7 



IGca^MiSMpi^aSHaS C2^M^'^Mp^o?a^ 2c2^M^'^Mp^a^a^ 
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4 // 2 /6 o 3 / /7 ryn 4 2 /8 



+ 



4c23Mi2Mpi^a4^2a6 2c2M^Mpi^a^n^a^ QAc2^M^'^Mpx^a^U'^a^ 

nc^^Pr,?^'^ ScsV^n^TT^^^^ 7p„,2a2 



15c3^PmTT'^f 5WpmTr'^ ISSCg^p^Tr'V 



4c2M9MpiVHa6 C2M6Mpi''aa4 4c2M6Mpi^aa4 

15c3V^7r^V^ 3c3^aV^7r^" 3c3^Pn.a^V^" 

8c2M6Mpi^aa4 4c2M6Mpi'^a!2^2a4 2c2Mf^Mp{^a^H^a'>^ 

csVm^a'vr'^ cs^Pm^aV'^ IScsp^vr'^/ 



2c22M9Mpi^a4H2a4 8c22M9Mpi^a3H2a4 2M3Mpi4aHa2 



4M6Mpi^^2q4 4M6Mpi^q^2q4 8c2^M^'^Mpi^a^n^a^ 
2703"^ Pm'^''^'H' ^ blcs^pma'Tr''^ IbCsPmT^'^f 



Wpm^Tl''' ^ 29C3Vm^7r'^ iWpma'TT'^f 



+ 



2c2^M^Mpx^a^Ha'^ 1^02'^ M^i^ a'^Ha'^ Ac2M^Mpi^a'^'Ha'^ 

CS^Pm^TT'^f 9C3V^V^V 3c3Pn^«V^^ 

2c2'^M^Mpi^a^na'^ %C2'^ M^i^ ama"^ M^Mpi^a'^V?a^ 



2M3Mpi4aH2a2 %C2M^Mp^am^o? 2,2c2M^Mpfa'^V?a? M^MpCUa? 
2c2^M^^Mpi^a^H'^a^ ?,2c2M^Mpi^aV?o? M^Mpi'^aHa^ ^ ' ' 



4Mpi2 128c23Mi8MpiiOa3^2aio 16c2^M^^Mpi^a^n^a^ 

33C35VV^2 3C35^^V>nV 15C3^7r'^p„,^ 

32c2^M'^^Mpi^a^'Ha^ 16c2^M^^Mpi^a^'Ha^ Uc2^M^'^Mpi^a^V?'a^ 

Scs^TT'^Pm' W^'^Pn^f W^'^Pm''^' 

8c2'^M^^Mpi^am^a^ 8c2M^Mpi^aa^ 8c2'^M^Mpi^a^'H^a'^ 

33C3^TT'^Pm^ Wi^'^pjf iW^'^pJ 

lQc2^M^Mpx^a'^na* 8c2^ Mpi^ a'^W 32c2M^Mpx^a'^rPa? 

3c3^7r^Vm^ 3pm^a'^ 
32c2M^Mpi^aV?o? SMpi'^H^ 



(C.4) 



^ 69C27r^^Pr»2 C27r^^p^2 ^ 19c27r^^pnt2 ^ 9C2«V^p,n2 

^^^^'^^ SMpi^^^ 8Mpi6a2^4 4Mpi6qW4 2Mpi^'H'^ 
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Mpi^n^ AMp\'^V? 2Mpi*n^ Ulpi^n^ 16c2^M^*Mpi^^a^n^a^^ 
8c2^M^^Mpi^'^a^n*a^^ 32c2^M-^^Mpi^^a*n^a^^ l&C2^M^^Mp]^^am'^a^'^ 

1227c3^y^V^2 165C3V^Vn^V 62303^^^^^^ 

32c23M2iMpii2^3^3ai4 + 32c23M2iMpii2^3^3ai4 + QAc^'^ M^^ M^i^"^ a^"^ o^'^ 

o 6 /12 // 2 o 7 /13 3 / loo 6 /12 2 /2 

_ 3c3°7r' q; Pm 3c3^7r' Pm « 183c3°7r' « 

2c23Mi8Mpii0a4^4„i2 + 8c24M2iMpii2a5H4ai2 ^ 32c23Mi8Mp/0a5H%i2 

16c23Mi8Mp/0a3H4ai2 16c24M2iMpii2a4^3ai2 8c23Mi8Mp/0a4H3ai2 

135c3^7r^^V^2^7 5541c367r^^V^2 

8c24M2lMpil2a47{3al2 16c23Ml8Mpil0Q4^3al2 + 32c23Ml8Mpil0a37^2Ql2 

465C36^-^V„.V 87c36V^V,„V2 cs^tt'^^pJ 

8c23Mi8Mpii0a3^2ai2 + 32c23Mi8Mpii0a37^2Qi2 2c23Mi8Mpii2a'i^4(^io 

29c3'7r^"p„,^a^ 81c3^7r^^V„,2^/ 

16c23Mi8Mpii2a3H4o,io 2c22Mi5MpiiOa4^4aio + 8c22Mi5MpiiOa3^%io 

o 5 /9 // 2 / o 5 /9 2 /3 o'y 5 /9 2 /i // 

3c3°7r q; prn « 3c3'^7r pm a 27c3°7r Pm a n. 

2c^M^Hdp^(^n^ c^M^^Mp^ofiH^ 4c2=^Mi5Mpi^a4^4g 

1375c357r^"p^2 1107c35y^V^2 27c357r^'aV^2 

16c22Mi5Mpii°a3^3aio 16c22Mi5Mpii°a2^3aio 4c23Mi5Mpi«a4H3ai 

3c3^7r^^°p^3a^ 243c3^VV^^a^^ SOcs^tt'^ Pm^H' 

8c24Mi8Mp/°a5^3aio 8c23Mi5Mpi«a5^3aio + c23Ml5Mpi«a3^3(jlO 

29C3VV^V 261c3^7r^^V^V 3cs'7r'\" p^' f 

2c22Ml5Afpil0Q:3^3alO 16c22Ml5Mpil0a2^3(il0 2c23Ml5Mpi8a4^3(jl0 

6c35VV.„^a^V 27c3^VV„,^7^7 15c367r^^°p^3 

C23Mi5Mpi8a5^3aio 4c23Mi5Mpi8a3^3aio + 8c24M18Mpi10q4^2„io 

_^14c3VVnV 3c367r^^°p^3y IQScg^TT^ V„.=^a7 

C23Mi5Mpi«a4^2aio 8c2'iMi8MpiiOa4^2Qio + 4c23Mi5Mpi8a4?^2Qio 

1335c3^7r^V^^ 483c357r^V„.V UIcs'tt'^ p„,^ f 

8c23Mi5Mpi'^Q3^aio 4c23Mi5MpiSa3^oio 8c23Mi5MpiSa3^(jio 



+ 



375c3^7r^V^^ 120c3^7r^'p^V 33c3V^V.„V^ 

2c23Mi2Mpi6a3a8 C2^M^^Mpi^a^a^ 2c2^ M^^ M/' 

175c3^7r^^°p^2 601c3^7r^^°p^2 39c3V^^V^^ 

16c2Mi2Mpii0a3H4a8 32c2Mi2Mpii0a2^4a8 8c2Mi2Mpii0aH4a8 

Icii'- (\ pi„- 9(::i'- (\ p,„- 3c;i '- p„r(\ J- 

C22Mi2Mp,«a4H4a8 4c22Mi2Mpi«a37^%8 C23Mi5Mpi«a4H2aio 

n 5 /9 3 / lop- 4 /8 2 /2 o'v 4 /8 2 /2 

9c3''7r' pm a 135c3^7r prria 87c3^7r p^ a 

8c23Mi5Mpii°a4-^4a8 l6c22Mi2Mpi^a5H4a8 16c22Mi2MpiSa4-^4a8 
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4c23Mi5Afpii0a3-H%8 ^02'^ M^"^ M^f a^W C2^M^'^Mp\^a^n^a^ 



2c23Mi2Mpi6a5^%8 C23Mi2M/a6^3a8 c2^M^'^M^x^a'^n^a^ 



2c23Mi5Mp/°a3^3a8 4c22Mi2Mpi«a3^3(j8 2c23Mi2Mp,6a4H%8 

2095c3VVn' ISQcaV^V^^ IScaV^'aVn,^ 

16c22Mi2Mpi^a3^2„8 16c22Mi2Mpi«a2^2„8 2c23Mi2Mpi6a4H2a8 

30c3V^V„,^«^' 225c3V^'/o>„^7^^ 103c3V^V^^a7 

C23Ml2Mpi*^a5?^2a8 + 4c23Ml2Mpi6a3?/2a8 8C22Ml2Mpi8Q4?/3a8 
495C347r'Vm'/ , 615c3V'VmV 3C3Va"p„2^ 



+ 



8c22Mi2Mpi8a3^2a8 8c22Mi2Mpi8a2^2a8 2c23Mi2Mpi6a4^2a, 



C23Ml2Mpi*^a5^2a8 c23Mi2Mpi6a37/2a8 16c22Mi2MpiSa37^2a8 



117c3V>™^/^ ^ 3C3V'V™^1^7' 225c3V'V™2a' 



16c22Mi2Mpi^a2H2a8 4c23Mi2Afpi6a3-^2Q8 2c23Mi2Mpi6a4Ha8 
^ 5C347r'Vm=^ 1503^'^^=^ 63c337r'Vm'a' 



4c22Mi2Mpii°a3^4a6 8c22Mi2Mpii°a2^4a6 8c2M9Mpi«a4H4a6 
93c3^7r^V^^a7 ^ 77c33VV^^a^ 15c337r'Vm'a' 



2c23Mi2MpiV4Ha8 8c2M9Mpi**a3H4a6 8c2M9Mpi**a2H4a6 

3c3V'Vm'a7' 3c3%'Vm'«'H' ^ 6C3=^7r'Vm'«'H' 



C23Mi2Mpi6a4Ha8 2c22M9Mpi6a4-^4Q6 02'^ M^M^i^am^a^ 

2QWt:'^pJ 22c3^TT'^Pm^ 291C337r'Vm^ 



3c337r'^a'Vm' 3c347r'^p„3„/ ^ 6c337r'V„2^/2 



2c22M9Mpi^a4^3a6 4c23Ml2Mpi^a4^3a6 C22M9Mpi^a5^3(j6 
45C337r'Vm^a'^ 15C3=^7r'V™2^/ ^ g^g3^/5^//^^^2 



'ic2^M^Mpfa^n^a^ C2^M^Mpi^am^a^ C2^ Mpi"" a^V? 

C22M9Mpi6a2^3„6 8c2M9Mpi8a3^3„6 4c2M9Mpi«a2:^3a6 
51C3VV^V 3c337r'^a"pmV 3c337r'^p^2^/2^ 



8c2M^Mpi^aH^a<^ 2c2'^M^MpX^am^a^ C2^M^Mpi^a'^H^a<^ 



2c22M9Mpi6a3H3a6 2c22M9Mpi6a2^3(j6 4c23Mi2Mpi^a3H2a6 

81c3-V^V^'a^ _ 135c3^7r^Vn^'«^ Scs^n'"^ Pm"^ f 



2c22M9MpiS4-^2a6 2c22M9Mpi^a3H2a6 4c23Mi2Mpi«a3H2a6 
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^Wt^'%J c^'k'^'oI'pJ 



2M^Mp^o?n'^of leM^Mpi^aH^a^ 2c2M^Mj'a^n'^af^ 

WT^'^a"pm^ Wi^'^Pm^'a' 9c3^7r"pm^a'f 

2c2M^Mpx^am^a'^ Ac2^M^Mpi^am'^a'^ C2^M^Mpi^a^H^a<^ 

^ Ycs^TT^p^^a" Sc^^TT'^pm'a" 21C3^'V™V 



8c2M6Mpi^a4^%4 c2M^Mpi^a^n^a'^ 2M^Mp{'aWa'^ 

scs^'^p^^n' ncs^TT'^p^^n' , c^tt'^pJ/ 



4c2M6Mpi6a4^3Q4 2c2M6Mpi6a3H3a4 4.M^Mpi^H^o? 



2c2M^MpY^am^a^ 2c2^M^Mp^^am^a'^ C2M^Mpi^a^'H^a^ 



^C2M^Mp{'o?-W'af ^C2M^Mp{'o?U'^af %C2M^Mp{'a^V?a'^ 



8c2M6Mpi''a7^2a4 4c2M6Mpi''a3^2Q4 ^C2M^ Mpx^ g?V? a'' 

60c,^7T''pm'f 2W^''pm^f 2WT^"pJf 



C2M6Mpi^a^2^4 8c2M6Mpi''a2^2^4 8c2M^Mpi^an^c 
Ci^Tr'^Pm^ Scs'^Tr'^^pm^ Scst:'^ prn^ a' 



2c2M6Mpi^a2-H4a2 4c2M^Mpi^an'^a'^ AM^Mpi^a^'H'^a^ 
23c37r''Vm2Q,/ 3c37r'^pm2Q/ ^ 171c37r^^/)m^ 3pm^a2 



4M3Mpi^a2^4a2 M^Mpi^aU'^a? AM^Mpi^W'a? 2Mp{'n'^ 

^Csn'^Prri^ 27c3Tr'^Pm'^ 215c37r^^Prn^ 

4M3Mpi6a3H3a2 4M3Mpi6a2^3a2 4.M^Mpi^aH^a^ ' 
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